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Abstract 

A variational perturbation theory based on the functional integral approach 
is formulated for many-particle systems. Using the variational action ob- 
tained through Jensen- Peierls' inequality, a perturbative expansion scheme 
for the thermodynamic potential is established. A modified Wick's theorem 
is obtained for the variational perturbation expansions. This theorem allows 
one to carry out systematic calculations of higher order terms without wor- 
rying about the double counting problem. A model numerical calculation 
was carried out on a nucleon gas system interacting through the Yukawa-type 

potential to test the efficiency of the present method. 
PACS numbers: 24.10.Cn, 71.10.Ca, 03.65.Ca 



1 



Typeset using REVT^ 



I. INTRODUCTION 



The perturbative approach has been proven widely successful when the interactions be- 
tween the particles were weak[l-3]. However, when the correlations between the particles 
are strong, the perturbation approaches are less than satisfactory. In order to overcome such 
difficulties, many nonperturbative schemes including variational calculation have been devel- 
oped [4-6] . Although the variational procedure provides a simple and convenient scheme to 
calculate the physical quantities, it has a serious drawback. Since, it is a basically one-shot 
process, systematic improvement over the obtained results is not easy to implement. 

In this paper, we present a functional integral formulation of variational perturbation 
scheme for many-particle systems. In fact, this concept has already been used in connec- 
tion with various research interests such as lattice dynamics, relativistic field theories, and 
quantum mechanics [7-12]. However, to the authors' knowledge, it has never been applied 
to nonrelativistic many-particle systems. For this purpose, we first derive the vatiation prin- 
ciple for the partition function and the thermodynamic potential in the functional integral 
formalism using Jensen-Peierl's inequality [13,14]. Then, this principle is used to obtain a 
variational action for the functional integral. The difference between the true and the vari- 
ational action provides a renormalized perturbation. Application of the Wick's theorem on 
this renormalized perturbation leads into cancellation of intracell joining diagrams. It will 
be shown that this modified Wick's theorem allows us to carry out systematic improvements 
on the variational result without worrying about the problem of double counting. As a test 
on the efficiency of the present scheme, we calculate up to the second order contributions to 
the ground state energy of interacting nucleon gas through the Yukawa-type potential using 
both the conventional and the present variational perturbation scheme. 
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II. THE VARIATIONAL PRINCIPLE IN THE FUNCTIONAL INTEGRAL 

APPROACH 



The variational principle for the partition function and the thermodynamic potential can 
be obtained using Jensen- Peier Is' inequality [13,14], which is given as follows for Hermitian 
operators Aq and Ai ; 

^j.g-(^0+^l) 



Ao - 



where 



Tre-^o 



> e-<^i^° , (1) 



Tre-^o 

When ^0 ^-nd Ai are normal-ordered operators composed of creation and annihilation op- 
erators, Aq = : Ao({aJj, tta}) '■ and Ai = : Ai{{al^, a^}) : , the trace operation is expressed in 
the coherent state functional integral form as follows [3], 

Defining the generating function K[f],r];AQ\ with source variables fj and r], 

K[fj,r]; Aq] = j D[ijjilj]e~ •i' '^'^^^'^'^"^'^^^^^"^'^^'^^^^^'^"^'^^'^^ , (3) 

we can express the above traces as follows, 

Tre-^o = /r[7y,r?;/lo]|^,,=o , (4) 
Tre-(^o+^i) = e-^'^^^'^^^^'^y^K[fj,r];Ao]\fi,r,=o , (5) 



(6) 



where ( is — 1(-|-1) and r] and r] are Grassmann (complex) variables for fermion (boson) field. 
Using the shorthand notation. 



dVa{r) ' dric^ir) J / /^^ K[r), rj; ^0]^^,^= 



(7) 



we can express Jensen- Peier Is' inequality in the functional integral expression, 



^0 ~ 



-(/'^-^i({<a;^'awfFj}))Ao 



(8) 



In order to apply this inequality to the partition function, we rewrite the action S in the 
form S — So{X) + {S — -S'o(A)), where So{X) is a variational action and A is a variational 
parameter to be determined through the variational calculation. Then the partition function 
is given by 



/ 



D 



= e 



D['iljilj]e 

(5[C 



-So[^,i>;\]-{S[-,PM-So[^,i':X\) 



d-q ' df} 1 



> e 



The thermodynamic potential can be expressed likewise 



(9) 



Q < -^In / D 



+ 



1 



s 



C 



d_ d_ 
dr]^ dfj 



or) or) 



(10) 



So 



We can obtain the variational result by maximizing the fourth line of Eq.(9) or minimiz- 
ing the righthand side of Eq.(lO) with A. Using the above expression, we now obtain the 
expressions for the variational perturbation process for the partition function and the ther- 
modynamic potential as follows, 

[^^] g-^o [V'.V': A] g- ( ('S'[C^ . ^] -'S'o [c^ , ^ ;a] ) ) 
g "( ^ ' ^ ]" ["^ ^ '^l ) + ( ^ '■!;]" ["^ -^^l ) s„ ^ 



So 



(11) 



Q=-lln/D 

-4ln/^-(^[^ 



dr) ' dfj 



'c^ ^-A 
or) or) 



J / So 



A. AA 

dri' drj] 



(12) 



P \ I So 

The second lines of Eq.(ll) and (12) are the perturbations based on the variational basis 
and, thus, provide systematic improvements over the variational result. 



III. VARIATIONAL PERTURBATION ON INTERACTING FERMION 
SYSTEMS AND MODIFIED WICK'S THEOREM 



The above variational perturbation result on the thermodynamic potential can be applied 
to any strongly correlated systems. In general, the variational result corresponds to the mean 
field value[15,16]. With the expressions given in Eq.(ll) and (12), we can readily carry out 
systematic improvements over the mean field result, since the renormalized perturbation 
term is generally much smaller than the original one. To see this point clearly, we apply 
the principle to an interacting fermion gas system through two-body interactions. In the 
process, we derive a modified Wick's theorem, which removes the problem of double counting 
of higher order terms. We also obtain explicit expressions for the thermodynamic potential. 

The partition function for two-body interacting fermion systems is generally given by 



where the subscript stands for all possible quantum numbers and an imaginary time r ; 
i = (kar). The integration or summation symbols are omitted by the summation convention. 
G° is a bare Green's function matrix and V is an interaction tensor. For example, for free 
fermions, G° has a form; 




(13) 




,-1 



{kar\G" 'ik'a'r') = (9^ + /2m - ii)5{t - T')6kk'5, 



(14) 



The partition function can be rewritten using the variational Green's function G ; 




v,v=o 



(15) 



Here we used the Gaussian integral for the Grassmann variables [3]; 




(16) 



Using the shorthand notation of Eq.(7), the partition function is expressed 
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, / rz-iQ /-'-J^N a a ir a o o a \ 

2^ _ ginDetG / ^'^a '^ij ) dm dfjj ^ijkiBrn drij ank ani \ (^yj^ 
Using the variational principle discussed above, we obtain for the partition function 



Z > Z — e ^ »3 »3 '\9vi ^"^i I Q 9vk I Q [V^ 

It is now necessary to obtain the Green's function, G, which appears in the variational 
process. For this purpose, we consider the variational thermodynamic potential. 



^<^= - ^InDetG-i - -{G%" - Q-^) '"^^ 



3 I G 



1^ Id d d d \ 
p '^''^ \ dr]i dr]j dfjk dfji / 



G 



1, „ 1 



= --lnDetG-i--(G°: -Gt^)G,, 

+ -^Vijki{GkjGii — GkiGij) . (19) 

Minimizing Cl against the variational Green's function G, we obtain 

Gba - GC + mbaj - Vuaj)Gji = . (20) 

In deriving the above result, the symmetry property of Vijki — Vjuk is used. Defining 

(6|S|a) = 2{Vibaj - Vbiaj)Gji , we obtain 

C-i _ G"'' + E = . (21) 

Here, we note that G corresponds to the self-consistent Hartree-Fock Green's function for the 
interacting fermion system. Using this variational Green's function, we rewrite the partition 
function as follows, 

2 ^ ^lnDetG-i+(G°^. ' -G'.^X^ g|- glj^^-Vijfcj <^ ^ g|- 

/ (co'^ ^n-^f _S--d / _S__d_\ \ _v- , , f S- S- S- -d /_8__9__S__2_\ ^\ 

•(e ^ij >\avidfij \dvi9fij/^J "^^kiyarndr,. Bfikdrii \ 877, Srj^- 9*7fc 9*?! / q M ^22) 

The second line in the above equation is the renormalized perturbation contribution and 
can be simplified as follows, 
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ydrii drij \dr]idr]j j ^) 
d d d d I d d d d 



ijkl 



drji drjj drjk drji \ drji drjj dr]k drji j 



( d d \ 

— (Ykijl + yiklj — Vikjl — ykilj) ~ j ^''^ 

( d d d d \ 

—Vijkl — T. T— T— GkjGli + GkiGlj 

\dr]; dr]j dr]k dr]i ' J 

\Or]i dr]j dr]k or]i or]j dr]k or]idr]i 

d d d d \ 

+Gij— — -—z- + Gki-^ — + GuGkj — GijGki 
■'drjidrjk drjj drji J 

^-Vrki(—— — —]' (23) 

'^^"^ \dr]i dr]j dr]k drji j 

In deriving the above relation, we used the symmetry of the potential Vijki and Eq.(20). 

Here the primed operation (• • •)' is defined to be 

fdddd\'_/dddd ^38 ^38 

\ 8r]i 8r]j 8fik 8fii J ~ \ 8r]i 8r]j 8fjk 8f]i 8r]j 8fjk 8r)i 8fji 

d d 8 8 \ 

TTz- + Gfcj- — + GliGkj — GijGki ] ■ (24) 
'8r]i8r]k 8r]j 8r]i J 

Using this primed operation and Eq. (21), the partition function is greatly simplified. 

1 / -V fS-S. a d \'\ 

1 ^ g-trlnG+itrEG / g "^H^iydrn dvj dfik 9m J \ ^25) 

In order to carry out the perturbation, it is necessary to evaluate terms like 

/f_d 8 8 d_\' /_9 8 8 d_\'\ 

\ \8rii, 8r]j, 8fik, 8fji, J 8fik^ 8fii^ ) j ^ 



8 8 8 8 \' ( 8 8 8 8 \' 



ri=ri=0 • 



(26) 



\dVii drjj^ drjk, 8rii, J \8rii^ 8r]j^ 8r]k^ 8r]i^ J 
With the prime defined in Eq.(24), we cannot apply Wick's theorem directly to the above 
expansion. Therefore, we first study contraction rules under the primed operation. We 
rewrite Eq.(24) as follows, 

( 4 4 4 4 ) ' = ( ^ ^' ^ n - ( )-{irki)-{iJ7i) 

I — I 

-{ijkl) + {ijkl) + {ijkl), (27) 
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where we used shorthand notations; ^ i and ^ — > i, and the connected pairs represent 
rn 

the contractions; i j = Gji. We call one (■ ■ ■)' a unit cell and mark a contraction with other 
unit cell elements with dots and a intracell contraction with lines. With these notations, 
terms in Eq.(27) appear 



ijkl) = {ijkl) + {ijkl) + {ijkl) + {ijkl) 



r~^_ I r~L L r~r~L L 

+ {ijkl) + {ijkl) + {ijkl), 



I ~"i I • ^ L I rn i 

.(^i j k I ) ^ -{i j k I ) - {i j k I ), 



-{i j k I ) = -{i j k I ) - {i j k I ), 



r—\_ . rrr i—t-L L 

-{i j k I ) = -{i j k I ) - {i j k I ), 



-( l)^-{iTk'l)-{ijk I). (28) 

Now adding all contributions within one cell, we find that all intracell contractions cancel 
and disappear, thus making the perturbation with the variational basis simple and straight- 
forward. The final expression for each primed operation in Eq.(26) appears 

(^^^A.)' ^ ( i j ^l) (29) 

V drji dT]j dr)k dr)i J ^ ■' ' ' V y 

Therefore, the contraction rule for higher terms now becomes 

( _d a a d_\' . . . /_d d d 9_y^ 

\drii^ dVii ^^h ) \9Vi„ drij„ dfjk„ dm„ J / 

'{ ii ji kih) - ■ ■ {in jn K In )] ^ 

(30) 



G 

Sum of all contractions between different cells. 



This modification of Wick's theorem represents nothing but a requirement that the selfenergy 
contributions which are included in the variational Green's function should not contribute 



again in the perturbation calculation. This new rule allows one to carry out higher or- 
der calculations without worrying about the double counting problem. The perturbation 
calculations are easily tractable by the diagramatic method. The unit cell with an inter- 
action, Vijki (^^^^) is shown in Fig.l (a). The modified Wick's theorem stipulates 
that diagrams with intracell contractions do not contribute to the thermodynamic potential. 
Examples of such diagrams are shown in Fig. 1 (b) . 

Since the linked cluster theorem is equally applicable in the present formulation, the 
thermodynamic potential can now readily be written 

n = ^trlnC - ^trEG - V e'^'^-'i^i^^^) 

where C indicates the connected diagrams only. The first two terms are the variational result 
and the last is the renormalized perturbation. For example, the second order perturbation 
terms are given as follows; 

11 I ( d d d d\' ( d d d d \'\ 

— —-^VijkiVi'j'k'i'Gi'iGk'jGkj'Gii' 
+iv^,.,V.,,„,G,.G.,G..,G,, , (32) 

which correspond to the diagrams in Fig. 2 (a). 




IV. A MODEL CALCULATION 

A sample calculation was carried out to test the efficiency of the present method us- 
ing a three dimensional nucleon gas system interacting through the two-body Yukawa-type 
potential. The model Hamiltonian is given by 

H^Yl efe«Lafe^ ^\ ^ 'v{q)al^al^,ak'-qa'ak+qa , (33) 

kcr kk'qaa' 

where 
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Here, m„ is the nucleon mass and 7 and A are model parameters which should be fixed for 
numerical calculations. The primed summation indicates that the q — Q term is renormalized 
away because this divergent energy term has the same contribution for both the conventional 
and the variational perturbation methods. With this model Hamiltonian, we calculate the 
ground state energy at zero temperature using both the variational perturbation theory and 
the conventional perturbation theory up to the second order. The bare Green's function and 
the interaction of this model is expressed as 

G--" ^ GCk'a'r' = {dr + 6. - ' r)6,k'Saa' , (35) 

^ijkl ^ ^'^l<Tiri,A;2<T2T2,A;3(T3T3,fe4Cr4T4 

= ^V{\k4 - ki\)Ski+k2,k3+kJ'Ti^JcT2a3S{Tl - T2)S{t3 - T4)5(t2 - Ts) , (36) 

where n is the chemical potential. Denoting GkarMr' Gkai^ — t') for a homogeneous 
system, the equation of G, Eq.(21) is expressed as follows, 

GUr r) = Gl{r -r')- f dT"Gl,{T - T")J2v{q)G,^,,{0)GUr" - r') . (37) 

Fourier transforming of this equation with respect to time, we easily obtain the variational 
Green's function Gkaii^n) as 

Gkcr{i^n) = -. 7^ T , (38) 

where Un is the fermion Matsubara frequency, (^"'+-^)^ and ika is the renormalized nucleon en- 
ergy value given by — v{q)nk+qa- nka is the Fermi distribution function, ^0{if,}-f,)^-^ with 
the renormalized energy. Therefore, the variational result of the thermodynamic potential, 

^min is 



1 1 

^mm — — trlnG 7:trSG 

2j3 



-\ E In (l + e-^(^'=-S'<"(^)'^'=+--'^)) + V v{q)nk+,ank. . (39) 



^ ka ^ kqa 
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The first order contribution is absent due to the modified Wick's theorem and the relevant 
second order diagrams are shown in Fig. 2 (a), while conventional diagrams of the same order 
are also shown in Fig. 2 (b) for the comparison, where double lines represent the variational 
G and single fines the bare G°. The details of the calculations are exactly the same for both 
cases except that the double line Green's function has the one particle energy, efc^- instead 
of the bare energy, e^. For the nucleon gas interacting with the potential given in Eq.(34), 
at zero temperature, the energy ika can be evaluated analytically to be given by 



~ - T 
- 2r 



(40) 



I J 7rm„ V4/ TT V 8 ; V^f' (97r/8)V3^ 

where r is the average distance between nucleons which is defined hy V — ^nr^N with the 
volume V and the number of nucleons N of the system, kp is the Fermi wave vector and 
F{x, b) is given by 



Fix, = - H m-- ) — 



_i 1 — x t1 + x 

tan — ; h tan — ; — 



(41) 



Using the above prescriptions and following standard text book steps [1] , we calculated the 
total energy at zero temperature numerically. The result is shown in Fig. 3 and compared to 
the conventional calculation. Since the exact solution is not available, it is not possible to 
access the accuracy of the present method directly. It is known that, for the homogeneous 
fermion gas, the variational Hartree-Fock result corresponds to the standard first order 
calculation [18]. Therefore, a convergence test can be carried out by comparing the second 
order contributions from both methods with the Hartrcc-Fock result. The result in Fig. 3 
clearly shows that the magnitude of the second order contribution is substantially reduced 
in the variational perturbation calculation. For example, at r = 2.5 Fermi, the renormalized 
contribution is 60.2% of the bare contribution. This renormalization clearly originates from 
the higher order terms included in the variational perturbation process [12]. 

The simple Yukawa-type potential considered here allows an analytic expression for the 
variational one-particle energy ika- However, in real systems, the interaction between nu- 
cleons is more complicated and, thus, generally does not yield a closed expression. We note 
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that, although an analytic expression for Ike is convenient to carry out the calculation, a 
numerical form of e^o- does not pose a fundamental difficulty, once whole spectrum of the 
energy values is known. 

V. SUMMARY 

In conclusion, we have presented a variational perturbation scheme for many particle 
systems using the functional integral formulation. It is shown that the Green's function 
obtained through the variational process corresponds to the Hartree-Fock Green's function. 
A modified Wick's theorem for renormalized perturbation based on the variational Green's 
function is obtained. A model calculation with a nucleon gas interacting with the Yukawa 
type potential shows that the present formalism provides an efficient and convenient formal- 
ism to study strongly correlated systems. 

ACKNOWLEDGMENTS 

This work was partly supported by the Korea Research Foundation(99-005-D00011) and 
by the Brain Korea 21 Project. 



12 



REFERENCES 

[1] A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems (McGraw- 
Hill, 1971). 

[2] G. D. Mahan, Many-Particle Physics (Plenum Press, 1981). 

[3] J. W. Negele and H. Orland, Quantum Many-Particle Systems (Addison- Wesley, 1988). 

[4] E. Fradkin, Field Theories of Condensed Matter Systems (Addison- Wesley, 1991). 

[5] J. Frohlich, Non-Perturhative quantum field theory (World Scientific, 1992). 

[6] R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

[7] T. R. Koehler, Phys. Rev. 165, 942 (1968). 

[8] P. Cea and L. Tedesco, Phys. Rev. D55, 4967 (1997). 

[9] G. H. Lee and J. H. Yee, Phys. Rev. D56, 6573 (1997). 

[10] H. Kleinert, Phys. Lett. A173, 332 (1993); W. Janke and H. Kleinert, ihid., 199, 287 

(1995) ; H. Kleinert and W. Janke, ihid., 206, 283 (1995); Annals Phys. 266, 135 (1998). 

[11] J. Krzyweck, Phys. Rev. A56, 4410 (1997). 

[12] S. K. You, K. J. Jeon, C. K. Kim, and K. Nahm, Eur. J. Phys. 19, 179 (1998). 

[13] L. S. Schulman, Techniques and Applications of Path Integration (John Wiley and Sons, 
1981). 

[14] D. Ruelle, Statistical Mechanics - Rigorous Results (W. A. Benjamin, 1969). 

[15] H. S. Noh, C. K. Kim, and K. Nahm, Phys. Lett. A204, 162 (1995); ihid., A210, 317 

(1996) . 

[16] H. S. Noh, S. K. You, and C. K. Kim, Int. J. Mod. Phys. Bll, 1829 (1997). 
[17] P. Ring and P. Schuck, The Nuclear Many-Body Problem (Springer, 1980). 

13 



[18] D. Pines and P. Nozieres, The Theory of Quantum Liquids (Addison- Wesley, 1989). 



14 



FIGURES 

FIG. 1. (a) Unit cell to construct Feynmann diagrams, (b) Diagrams which do not appear in 
the variational perturbation calculation. The double line indicates the variational Green's function. 

FIG. 2. Second order diagrams (a) for the present variational perturbation scheme and (b) for 
the conventional perturbation theory. Note that the last diagram in (b) does not appear in (a) due 
to the modified Wick's theorem. 

FIG. 3. The ground state energy per particle calculated up to the second order; (a) solid line 
represents the present variational perturbation result, (b) dotted line the conventional perturbation 
upto the second order, and (c) dashed line the Hartree-Fock result. Here, is the Compton 
wavelength of the pion [17]. 
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